Particle focusing in oscillating dissipative billiards 
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We develop and analyze a scheme to achieve both spatial and energetic focusing of an ensemble 
of neutral particles which is based on an oscillating billiard with frictional forces. The interplay 
of two competing mechanisms, acceleration due to collisions with the oscillating billiard walls and 
decceleration caused by friction, leads to the emergence of attractors in phase space. Their specific 
properties, i.e. spatial localization and energy spread, can be controlled and tuned by varying e.g. 
the frequency of the time-dependent billiard. 
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Spatial and energetic focusing of neutral particles is 
an indispensable prerequisite for many experimental se- 
tups such as atomic [l| and crossed molecular beams [2^ , 
microparticles in a 'lab- on- a- chip' [sl or in ultracold neu- 
tron precision measurements y. Typical devices used 
in order to spatially localize particles are lenses, mirrors, 
nozzles or collimators. The desired energies are usually 
achieved by means of velocity selectors or, in the case of 
neutrons, by moderating media. In this work, we propose 
a novel scheme for spatial and energetic focusing of an en- 
semble of neutral particles by using a remarkably simple 
billiard setup. Necessary ingredients are oscillating walls 
and frictional forces acting on particles in between colli- 
sions with the walls. Their interplay leads to a state of 
dynamical equilibrium, i.e. periodic and quasi-periodic 
attractors in phase space, whose energy and spatial lo- 
calization can be tuned by varying e.g. the oscillation 
frequency of the billiard. 

It is well-known that time-dependent billiards [Hl-lsl can 
be used to study non-equilibrium phenomena like Fermi 
acceleration IqI, 1iO| and anomalous diffusion in momen- 
tum space |lll4l3|. Moreover, dissipation can be intro- 
duced to billiard systems either by a restitution coef- 
ficient accounting for inelastic collisions 14-16|, or by 
adding d rag forces between collisions (e.g. Stokes' fric- 



tion) llTHl 



Here, we investigate the particle dynamics in a 2D har- 
monically driven elliptical billiard with a frictional force 
between collisions. There are two competing mechanisms 
in the presence of driving and dissipation: Firstly, the 
driving accelerates the particles and secondly, the fric- 
tion slows them down, i.e. we expect a dynamical equi- 
librium where both processes balance each other. We will 
demonstrate that the emergence of attractors focuses the 
ensemble onto certain parts of phase space. The posi- 
tion of the attractors in phase space can be controlled 
by tuning the parameters of the driving. Consequently, 
the dissipative driven elliptical billiard can be used as a 
tool for particle focusing in space and energy without any 
imaging or focusing devices such as lenses, collimators or 



monochromators. 

The boundary Bit) of the driven elliptical billiard 
is given by B(t) = {(a(t) cos 0, 6(t) sin^) |0 < (/> < 27r}, 
where a{t) = ao + C sin(cjt) and h{t) = ^ C sin(cjt) 
are assumed to be harmonic functions, where C > 0, 
uo and ao,6o are the driving amplitude, the oscillation 
frequency and the equilibrium values of the semi-major 
and semi-minor axes, respectively. Due to friction, par- 
ticles in-between collisions with the moving boundary do 
not travel ballistically, but they still move on straight 
lines: v = kv v{t) = Voe~^^ and the position 
of a particle (in-between collisions) is simply x{t) = 
Xq -\- ^ [l — e~^^). We focus on the case (sufficiently 
small values of /c), where a stopping of the particles in- 
between collisions occurs hardly ever. 

The dynamics can be described by an implicit AD map- 
ping, using the variables (^n, (/^n, ci^n, '^n), where Vn = 
I'^^nl, £,n = cjtn(mod27r) is the phase of the boundary 
oscillation and is the angle between Vn and the tan- 
gent of the boundary at the n-th collision (see Fig. [1]). 
The time tn+i of the n + 1-st collision is determined by 
the smallest t^+i > tn that solves the implicit equation 
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where At = tn+i — tn and Vn = (v^^yy). The equations 
for 0n+i7 <^n+i and Vn-\-i are the same as for the elliptical 
billiard without dissipation, see Eqs. (9)-(12) of Ref. [l3| . 

Without dissipation, the driven elliptical billiard shows 
FA [lll-[l3|, i.e. the ensemble averaged modulus of the ve- 
locity {v){n) grows with increasing number of collisions 
unboundedly. With friction, this does not hold anymore. 
The evolution of {v){n) as a function of the number of 
collisions n is shown in Fig. [2l Initially (n < 10^), 
{v){n) increases, roughly according to a power law, just 
like in the elliptical billiard without dissipation. This 
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FIG. 1: (Color online) Boundary of the elliptical billiard in 
coordinate space at two different phases ^ together with a 
typical period four attractor (ao = 1.5, 60 = 1, C = 0.25, c<; = 
1). At the single crossing point, the two fluxes $1 and $3 
intersect. 



transient increase is followed by a plateau-like behav- 
ior (10^ < n < 10^) and a subsequent second increase 
(10"^ < n < 10^). Finally (n > 10^), the mean veloc- 
ity {v){n) saturates to a constant value, i.e. there is no 
FA in the long-term dynamics (this has been checked nu- 
merically up to n ~ 10^). Naturally, the saturation value 
depends monotonically on the dissipation parameter k 
and increases with decreasing k. 

Intuitively, the absence of FA can be understood in the 
following way: Between two successive collisions with the 
oscillating billiard boundary, a particle will slow down by 
a certain amount Avf due to the presence of friction. For 
large the velocity loss Avf becomes also large since it 
is proportional to v itself. On the other hand, a particle 
can gain momentum upon a collision with the boundary 
if the billiard is contracting. However, the velocity gain 
Avc < 2uoC upon a single collision is independent of v. 
Overall, this means the driven billiard cannot support 
velocities larger than a typical scale and consequently 
fast particles will slow down until they are in the velocity 
range of the attractors. 

Effectively, this constitutes high-lying energetically 
forbidden regions in phase space for particles starting at 
low velocities. In the beginning of the time-evolution, the 
particles spread diffusively in phase space, which leads to 
the transient (n < 10^) increase of {v){n). At the same 
time, the standard deviation cr^(n) of the velocity distri- 
bution pn{v) increases as well, see Fig. [2l The diffusion 
stops when the ensemble is uniformly spread over the 
accessible part of phase space, resulting in the plateau- 
like structure of {v){n) between 10^ and 10^ collisions. 
Concurrently, cr^(n) saturates. 

The reason for the second increase of {v){n) (10^ < 
n < 10^) is the following: Due to the presence of dissipa- 
tion, certain attractors emerge in phase space. These are 




FIG. 2: (Color online) {v){n) (dashed) and its standard devia- 
tion av{n) (solid) (1000 particles with vo = 0.1, and ^0, 0o, cxo 
distributed uniformly; C = 0.25, ao = 1.5, 60 = 1, /c = 10~^). 
The inset shows the distribution pn{v) at n — 10^ (dashed) 
and at n = 10^ (solid). 



either stable limit cycles, corresponding to quasi-periodic 
motion, or stable spiral points, corresponding to periodic 
motion. For the given set of parameters, the velocity of 
particles on (or very close to) such attractors is roughly 
between 0.5 and 1.5. Step by step, the particles are cap- 
tured by these attractors, hence the standard deviation 
ay{n) decreases in this (10^ < n < 10^) regime. Since 
the attractors lie rather close in momentum space to the 
maximum possible energy [I^, the capturing of the par- 
ticles onto the attractors is associated with an increase 
of {v){n). Eventually (at around n ^ 10^), almost all 
initial conditions are captured by these attractors and 
consequently {v){n) and cr^(n) approach constant values. 

To study the capturing process of the particles onto the 
attractors further, we investigate the (collision resolved) 
evolution of the phase space density 



— r dvp{^,(t),a,v,n)dv, 
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(2) 



where p{^, a, ^, n) = ^ E^L^ ^(^ - - ^m^^ - 

a2)S{v — vf) is the phase space density of an ensemble 
of particles at collision number n. 

In the beginning (Fig. [3^), the ensemble is distributed 
almost uniformly over the whole (j) x a plane. With in- 
creasing number of collisions (Fig. [Sb) some particles get 
focused in the immediate vicinity of the limit-cycle type 
attractors. These attractors correspond to particles skip- 
ping along the elliptical boundary almost tangentially, 
i.e. they are whispering gallery orbits and can be ob- 
served in Fig. [3t) as wavy hues at a ^ 0, tt. After a very 
large number of collisions, almost all particles are cap- 
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tured by attractors (Fig. [Sh), i.e. the phase space den- 
sity is zero outside the regions covered by the attractors. 
Now, both types of attractor, periodic orbits (bright dots 
(22! ) and hmit-cycles are populated. 

There co-exist several periodic attractors with differ- 
ent periods. These periods are typically between two and 
thirty (number of collisions with B{t)). Exemplarily, the 
trajectory of a period four orbit in coordinate space is 
shown in Fig. [TJ together with the elliptical boundary at 
two different phases ^. Since the orbit closes after four 
collisions, the net energy gained by the four boundary col- 
lisions has to exactly match the energy that is dissipated 
in between these collisions. Furthermore, the orbit has 
to close after four collisions in the other three variables 
(^, 0, a) as well. Fulfilling all these conditions simultane- 
ously is highly nontrivial, being the reason why there are 
usually only a few (up to ten) dominant periodic orbit 
type attractors for a given set of parameters. 

The presence of just a few dominant periodic attractors 
renders the driven, dissipative elliptical billiard an ideal 
system in terms of controllability. To remove the vast 
number of limit-cycle type attractors that undermine this 
controllability, it is sufficient to introduce a small hyper- 
bolic element (e.g. a small hump) to the boundary of the 
billiard at = and (j) = n. Particles on stable limit- 
cycles will hit this hyperbolic element and get scattered 
back into different parts of phase space. Eventually this 
process will drive all initial conditions onto the periodic 
attractors. Note that we choose ^ = 0, tt for the positions 
of the hyperbolic elements since generically the periodic 
attractors do not have collisions with the boundary at 
these values of (j) (see Fig.[3h). The resulting velocity dis- 
tribution pn{v) in the long-term asymptotics (n > 10^) 
is shown in the inset of Fig. O Specifically, there is one 
dominant period four attractor (see Fig. [T]). Each of the 
four branches of this attractor is associated with a dif- 
ferent velocity {vi = 1.30, 0.95, 1.33, 1.50, i = 1 . . . 4). As 
a consequence, pn^io^i"^) of the ensemble of particles is 
not smooth, in contrast to pn=io^{^)^ but shows distinct 
peaks at exactly these four values of Vi. 

By tuning the parameters of the system, especially the 
frequency a;, the amplitude C or the friction constant /c, 
the position of the attractors in momentum space, and 
thus p{v) can be controlled in a desired manner. Since 
the frictionless billiards exhibits Fermi acceleration, its 
dissipative counterpart can possess attractors at arbitrar- 
ily large energies. Furthermore, we can populate a cer- 
tain periodic attractor exclusively by adding additional 
small humps to the elliptical boundary at positions where 
particles on other (unwanted) periodic attractors hit the 
boundary. 

With the above described procedure, particles that are 
not on the attractor of choice are reinjected into different 
parts of phase space until finally all particles are captured 
by a desired attractor. Thus, ensembles of particles con- 
taining few distinct velocity components can be created. 
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FIG. 3: (Color online) Collision resolved evolution of the 
phase space density (logarithmic colormap) /)ni,n2 (0, <^): a) 
10^ < n< 2-10^, b) 4-10^ < n< 10^ and c) 2-10^ < n< 10^ 
The ensemble gets focused with increasing number of colli- 
sions more and more onto attractors (bright dots and wavy 
lines in panel c)). 

By temporarily providing a small opening of the billiard 
at the collision point of a certain branch of the attrac- 
tor with the boundary, a monoenergetic spatially focused 
ensemble can be coupled out of the billiard. 

Finally, we investigate the impact of particle-particle 
interaction on the focusing process onto the periodic at- 
tractors. If we assume a short range i.e. contact inter- 
action, there are now two competing processes. Firstly, 
'background particles' (particles that are not yet on this 
attractor) get captured onto the attractor with a certain 
rate Tc- Secondly, particles on the attractor get scattered 
out of this attractor by two different mechanisms: 

1. Particles interact with background particles, result- 
ing in a scattering rate r^. 

2. At these crossings points of the attractor in coor- 
dinate space, captured particles interact and thus 
scatter with particles on the same attractor with a 
certain rate r^. 

This can be modeled by the following system of coupled 
differential equations for the population Na of a certain 
attractor in the long-terms dynamics 

Na = TcN, - TsNl - nN.Na (3a) 
Nb = -rcNb + rsN^ + nN.Na, (3b) 
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where is the number of background particles. The 
capturing rate Tc, which of course depends on the spe- 
cific attr actor, is determined from our numerical simu- 
lations. In the following, we consider exemplarily the 
period four attractor shown in Fig. [T] employing small 
humps at (/) = 0,7r to remove the whispering gallery at- 
tractors. To calculate the scattering rate rg, we model the 
(short range) interaction by a typical scattering cross sec- 
tion a. At the crossing point in coordinate space, there 
is an incoming flux <l>i of particles along branch 1 and 
also a flux <l>3 along branch 3. The number of particles 
on branch 1 that scatter within a unit time interval at 
the crossing point is given by TVi = pi^i^ where pi is the 
probability for a particle on branch 1 to interact with a 
particle on branch 3. During the time interval in which a 
particle on branch 1 travels the distance a in order to pass 
the crossing point, a particle on branch 3 travels the dis- 
tance s = (TVs I '^1 • The scattering probability p\ is then 
Pi = (cr + sj^a/'^s. Since the fluxes obey <l>i = <l>3 = <l>, 
this yields N\ = ^l!^i±!^<|)2 ^-^^ the same expression 

holds for 7V3. The flux <l> can be written as <l> = Na/Ta^ 
where Ta is the time for one round trip on the attractor. 
Thus, the change of particles Na = —Ni — Ns on the 
attractor due to scattering at the crossing point and the 
scattering rate (compare with r.h.s. of Eq. ([3a|) ) are 
given by 



2(t{vi + vs) 
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From kinetic theory j20| one obtains for the scattering 
rate between background particles and particles on an 
attractor 



(5) 



where (vb) is the mean velocity of the background parti- 
cles and Aei = fr J^^ a(t)h(t)dt = 7r(ao6o + is the 
mean area of the elliptical billiard. 

The system of differential equations (|3|) possesses sta- 
ble solutions that approach asymptotically a constant 
value. Specifically, we obtain numerically rc = 3 • 10 ~^ 
and for (J = 2 • 10~^ the scattering rates are Vg = O.OSa 
and rb = 0.3a {{vb) = 1) yielding a high population of 
Na = 94.3%. As long as a is smaller than 10~^ao, the 
dynamics, especially the focusing of the ensemble onto 
attractors is not disturbed significantly, implying the pos- 
sibility of a controlled preparation of the ensemble. 

We demonstrated that combining dissipation and driv- 
ing in a billiard leads to the emergence of periodic attrac- 
tors in phase space onto which the originally uniformly 
distributed ensemble gathers. Our remarkably simple 
setup therefore allows spatial and energetic focusing of 
particles and might represent an experimental alternative 
to devices such as lenses, collimators and monochroma- 



tors. In principle, the billiard setup could be of macro- 
scopic or mesoscopic character. One possible realization 
is an atom-optical experiment, since there, the geometry 
of the billiard (a rapidly scanning blue-detuned laser) 
can be varied arbitrarily in time using acousto-optical 
modulators and dissipation can be introduced via opti- 
cal molasses (ilj . 
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